The Granular Character of Particle Rafts 
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We consider a monolayer of particles floating at a horizontal liquid-gas interface — a particle raft. 
Upon compressing the monolayer in a Langmuir trough, the particles at first pack but ultimately 
the monolayer buckles out of the plane. We measure the stress profile within the raft at the onset 
of buckling and show for the first time that such systems exhibit a Janssen effect: the stress decays 
exponentially away from the compressing barriers over a length scale that depends on the width 
of the trough. We find quantitative agreement between the rate of decay and the simple theory 
presented by Janssen and others. This demonstrates that floating particle rafts have a granular, 
as well as elastic, character, which is neglected by current models. Finally, we suggest that our 
experimental setup may be suitable for exploring granular effects in two dimensions without the 
complications of gravity and basal friction. 



A monolayer of densely packed particles floating at a 
liquid interface has many curious and potentially useful 
properties. At a coarse-grained scale, this particle-coated 
interface is reminiscent of an elastic sheet since compres- 
sion leads to a buckling instability 0, Q while in ten- 
sion it may also fracture[3|. Detailed experiments con- 
firm a transition from liquid-like to solid-like behaviour 
as the particle concentration increases Qj. At the scale 
of the individual particles, the complex shape of the 
liquid-gas interface has been shown to prolong the life- 
time of particle-coated bubbles[5|. Furthermore, col- 
loidal particles are well known to be capable of stabilising 
emulsions and enable the formation of dry capsules or 
colloidosomesQ. This rich variety of properties ensures 
that particle-coated interfaces have found application in 
diverse arenas from drug delivery[7[ to waste disposal by 
insects @. 

At the same time, an intensive and largely independent 
research effort has focused on understanding the stress 
state within granular materials. One of the classic papers 
in this field is by Janssen0,[I3| who considered the stress 
distribution in a granular silo. Janssen showed that the 
vertical pressure distribution in a silo filled with grains 
is not hydrostatic. Instead, the pressure saturates at a 
constant value which scales as the weight of grains over 
a vertical distance comparable to the width of the silo. 

In this Letter we study the buckling instability of a 
particle raft, and address the question of what causes 
the raft to buckle. Previous experiments on very small 
(colloidal) particles and theoretical calculations have sug- 
gested that bucklin g oc curs when the effective surface 
energy becomes zero[llJ. This is also the case for surfac- 



tant or lipid monolayers [FJ, |13|. In these experiments it 
is usual to give results in terms of the surface pressure 
IT = 7 C — 7, in which 7 C is the surface tension coefficient 
of the clean interface (i.e. in the absence of particles) and 
7 is the effective surface tension of the 'contaminated' in- 



terface. Thus the orthodox view is that buckling occurs 
when IJ = 7 C . Experimentally, it is conventional to mea- 
sure the surface pressure using a Wilhelmy plate: a strip 
of material that is wetted by the underlying liquid and is 
inserted into the surface. By measuring the force exerted 
on the strip by the liquid interface, the surface pressure 
LT may be inferred. Recently, it was suggested that this 
method of measurement may give distorted results when 
the interface is no longer fluid-like [3], e.g. when a par- 
ticle raft is very densely packed. Indeed it is well known 
that for viscoelastic or elastic layers, the surface pres- 
sure can be anisotropic, and is best regarded as a surface 
stress Here, we present the first systematic test of 
the orthodox view that II = 7 C at the onset of buckling 
and show that, as Pocivavsek et al. intimated [F4(, the 
surface pressure is not uniform throughout the raft (al- 
beit for a different physical reason). We find that the 
pressure measured at buckling by a Wilhelmy plate de- 
pends on the width of the trough and the separation of 
the barriers at this point. This observation cannot be 
explained using current elastic models and we show that 
it is a consequence of the granular stress state in the par- 
ticle monolayer. Crucially, this demonstrates that such 
systems have granular characteristics that have not been 
considered previously. 

Our experimental setup is a monolayer of relatively 
large, athermal, particles[27| at an air-water interface in 
a commercial Langmuir trough. Two different troughs 
were used: Minitrough (KSV, Finland) and NIMA 600 
(NIMA, UK). Particles are added to the air-water inter- 
face in a petri dish of area « 60 cm 2 until the addition of 
further particles would lead to interfacial buckling. This 
gives a reproducible 'unit' of a large number of particles 
which are then transferred to the trough. The pliolitc 
particles are so hydrophobic that none sink or are lost 
during this transfer. Before each run of the experiment, 
the interface is agitated manually to ensure that any par- 
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FIG. 1: The exerimental setup, (a) We use a Langmuir trough 
with barriers attached to a computer-controlled stepper mo- 
tor. The width of the trough, w, is varied by moving the 
strips beneath the barriers. The computer is calibrated to 
infer the surface pressure If from the force on the Wilhelmy 
plate (red rectangle) . In the configuration shown in the main 
figure it measures ILi. To measure II_l, the plate is rotated by 
90° in the (x, t/)-plane, as shown in the inset. Dimensions of 
the trough itself are for the NIMA 600. (b) Images showing 
the center of raft at three stages of compression. From left 
to right: initial state, onset of non-zero pressure and final, 
highly buckled, state. The scale bar represents 1 cm. 



tide clumps are broken up and the particles distributed 
evenly throughout the trough. The Langmuir troughs 
have two motorised barriers that are computer controlled 
so that the monolayer may be compressed symmetrically 
at constant barrier speed (the speed of compression is 
100 /nm/s|28j). A Wilhelmy plate placed at the center 
of the trough is used to measure (via supplied computer 
software) the surface pressure II as the monolayer is com- 
pressed. Two different orientations of this plate were 
used, as shown in fig. [T^,. In each experiment we use the 
sensor in one of these two configurations. We denote the 
surface pressures measured with the Wilhelmy plate by 
II_l and IX 1 1 depending on whether the plate is perpen- 
dicular or parallel to the mobile barriers of the Langmuir 
trough. We are able to vary the width, w, of the Lang- 
muir trough available to the monolayer by using specially 
cut strips of Delrin, placed in the subphasc flush with the 
water surface, as shown in fig. QJi. 

This experimental system is different from that con- 
sidered by Aveyard et al. ll| and others in three impor- 



tant ways. Firstly, in our system there is no long range 
repulsion between the particles (in fact, they are large 
enough that the attractive interaction energy from cap- 
illary forces is greater than thermal energy [l^). The 
only repulsive interaction is the close-range steric repul- 
sion. Secondly, the particles in this study are not perfect 
spheres and instead have rough shapes. Thirdly, our ex- 
perimental system is unique because of its facility to vary 
the aspect ratio of the trough. 

We measure the surface pressure as a function of the 
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FIG. 2: Examples of experimentally measured surface pres- 
sure isotherms showing II as a function of the distance be- 
tween the barriers and the sensor, d. II increases rapidly as 
d decreases but then exhibits a kink (marked by x in b). 
This kink corresponds to the onset of buckling. Isotherms 
are shown for a range of trough widths, w, and for different 
numbers of particles spread on the surface, (a) With a fixed 
number of particles but decreasing trough width, the value of 
H || measured at the onset of buckling decreases, (b) For fixed 
trough width but increasing number of particles, the value of 
XT || at buckling decreases. 



barrier separation, 2d. This produces surface pressure 
isotherms such as those shown in fig. [5] From these 
isotherms, we see that the surface pressure rises as the 
distance between the barriers decreases (i.e. as the mono- 
layer is compressed) . After the onset of a non-zero pres- 
sure, the pressure increases rapidly with further compres- 
sion until at some critical barrier separation, 2d c , we ob- 
serve a kink (marked by x in fig. ^p) . This is the com- 
pression at which the monolayer first buckles (monolayer 
collapse) . It is unlikely that a simple analytical form can 
be used to describe the measured isotherms. We focus 
instead on the measured value of the pressure at buck- 
ling, H buck . From the range of isotherms shown in fig. [21 
we see that Y[ buck depends on the width of the trough 
and also on the barrier separation at collapse, 2d Cl which 
in turn depends on the number of particles spread on the 
interface [29(. 

Several authors have expressed concerns that the sur- 
face pressure measured by the Wilhelmy plate may not be 
an accurate representation of the effective surface tension 
of the interface in the vicinity of the plate LJ, [l7| • To test 
this directly, we took photographs of the meniscus shape 
around the Wilhelmy plate. For a given value of the 
surface pressure measured by the Wilhelmy plate, we de- 
termine the interface shape expected for a pure interface 
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FIG. 3: Snapshots of the meniscus formed around the Wil- 
helmy plate compared to the theoretically predicted meniscus 
shape. Predicted shapes (solid curves) are obtained by solv- 
ing the Laplace- Young equation pl| with the value of the 
surface tension corresponding to the measured surface pres- 
sure: 7 = 7 C — n. 

with this effective surface tension coefficient by solving 
the Laplace- Young equation with zero contact angle at 
the plate [l^. Such a comparison is shown in fig. [3] and 
shows excellent agreement between the theoretical and 
experimentally observed meniscus shapes. We therefore 
conclude that the measurements of the surface pressure 
by the Wilhelmy plate are consistent with the observed 
interface shape. 

Having seen that the aspect ratio of the trough at col- 
lapse, 2d c /w, influences the surface pressure at the onset 
of buckling H buck , we now present experimental results 
for a range of trough widths and number of particles at 
the surface. The pressure H buck is defined as the kink in 
the surface pressure isotherm. To automate the detec- 
tion of this point, we determine the point of inflection of 
the isotherm, marked (*) in fig. [^L, and the point where 
the derivative falls to small values, marked (x) in fig.[2jD. 
In fig. [4] we show the results of these experiments pre- 
senting the pressure at buckling H buck as a function of 
d c /w. We performed 163 independent experiments, but 
for clarity of presentation we have averaged the results 
into bins that are equally spaced in d c /w. These bins 
range from d c /w = 0.1 to d c /w = 3.5 — a span of one 
and a half decades. Error bars represent the standard 
deviation of data values within each bin. We observe 
good collapse of data from a wide range of experiments 
in different troughs. This master curve suggests an expo- 
nential decay in H buck with d c /w. The prefactor of this 
exponential depends on the algorithm used to determine 
H buck but the decay rate is insensitive to this. We there- 
fore focus on this decay rate as the most robust feature 
of the results presented in fig. 0J 

We use ideas from the theory of granular elasticity [l9| 
to explain the exponential decay presented in fig. [4] Un- 
der the assumption that the material is elastic, and that 
there is no strain in the y-direction (perpendicular to the 
direction of compression, see fig. [Th.) we find that 

cr y = va x (1) 

where v is the Poisson ratio [2(|. We note that |T|) is 
simply the relationship posited by Jansenn between the 
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FIG. 4: Experimental data showing the pressure measured 
at the onset of buckling, n 6 " cfc , as a function of the aspect 
ratio of the trough at this point, d c /w. Results are shown 
for (a) Ilx and (b) II||. In each panel there are two datasets, 
corresponding to two automated measures of Yl buck from sur- 
face pressure isotherms: the inflection point (open symbols) 
and the point where the derivative falls to small values (filled 
symbols). Fits to the exponential decay of eq. ((4} are shown 
with characteristic decay parameters A = 1.8 and 2.3 for Ilx 
and A = 1.6 for both datasets of ELi. The prefactor in the 
exponential decay is sensitive to the definition of the buckling 
point and so is not considered here. (A power law decay is 
not appropriate judging from the corresponding logarithmic 
plot and would not yield finite U buch in the limit d c /w — > 0.) 

vertical and horizontal pressure in a silo0, [l(| : here, 
the constant of proportionality K — v 20\. Furthermore, 
we note that this relationship should hold only close to 
buckling since prior to this particles can, and do, accom- 
modate deformation by moving in the y direction. 

In our two-dimensional situation, a x = IX n and 
(Xy = II_i_. We now use a well-known balance of forces 
argument [2ll] to determine ITn (ar) . (Though LT|| may in 
fact be a function of y, we consider only its x dependence 
here.) To determine the behavior of H || we consider an 
infinitesimal slice through the trough (taken at constant 
x) with width w. In equilibrium, the difference in pres- 
sure between the two faces of the slice must be balanced 
by frictional forces acting at the walls of the trough, de- 
noted a X y, he.: 



Assuming a Coulomb friction law and a friction coeffi- 
cient [i we have that a xy — /iLIx = /wILi. We may then 
integrate eq. ([2]), which has solution: 

H]i (a?) = ccexpf— 2fj,i/x/w), (3) 
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for some constant of integration a. 

We observe experimentally (see fig. [Tp) that buckling 
occurs first in the vicinity of, and parallel to, the com- 
pressing barriers. Under the assumption that this cor- 
responds to a critical local pressure ILi = 7 C then the 
pressures measured by the Wilhelmy plate should be 

n|| bucfc = 7c exp(-2/W c /w), n ± buck = iyn {l buck . (4) 

In this relationship the constants /i and v are the only 
unknowns. Two values of v have been proposed for the 
case of perfect discs: v = 1/3 [53| and v = l/V^tH]. We 
measured the friction coefficient fj, directly; the angle of 
friction between Pliolite and the Dclrin strips was mea- 
sured to be 28° so that ju f» 0.53. Using these values 
we expect log(njj ucfc /7c) = -(d c /w)/\ with A = 2.83 or 
A = 1.63 using v = 1/3 and l/\/3, respectively. These 
predicted lines have a slope of the same order of mag- 
nitude as that measured experimentally. This suggests 
that the pressure decrease observed at the center of the 
raft relative to the value near the compressing barriers 
is indeed induced by granular effects. We note that 
taking v = 1/V3 w 0.577 gives a significantly better 
account of the experimentally measured slope than the 
value v = 1/3. Motivated by eq. {I}, we calculate the 
ratio of U^ ck (d = 0) and Hjj ucfc (d = 0) to determine 
an independent estimate of the Poisson ratio, v. This 
ratio gives v = 0.63 ± 0.05. v may also be estimated 
from the measured values of A. This procedure gives 
v = 0.52 ± 0.09. These values are consistent with one 
another and with the value for rigid discs [H. 

Previously, particle rafts have been modelled as elas- 
tic sheets [l|, 0, li , 14 1 . In this work we have shown that 
these approaches are insufficient — even very recent elas- 
tic calculations |14| do not predict the exponential decay 
reported here. An alternative interpretation of the last 
three data points in figure |4Jd is that H buck plateaus at a 
constant value for large values of d c /w. However, even in 
this instance the granular character of the raft remains 
clear. The granular character of such rafts is therefore, in 
addition to their elastic character, a vital ingredient that 
must be considered when modelling them in the range 
of industrial^ and biological^ settings in which they 
arise. 

Recently, rafts of bubbles and millimetric beads 24 1 
have proven useful for understanding shear-banding and 
glassy dynamics. In the same way, the present system of 
rough, floating grains may prove useful as a model system 
for studying the influence of friction in two-dimensional 
granular materials. Furthermore, using our system to 
study the jamming transition would remove the possi- 
ble effects of basal friction [2^] and allow confirmation of 
recent theoretical predictions for the influence of inter- 
particle attraction[26| by exploiting the small capillary 
attraction In addition, the Wilhelmy plate provides 



a convenient method of probing the stress state within 
the medium. 
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